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There exists in nature many examples of systems presenting self-limiting behaviour: population 
dynamics, structure engineering, Townsend's electron breakdown, nuclear decay in radioactive equi- 
librium, histeresis process, meteorological models, etcetera. In this work we call attention to the 
advantages the use of a variational formulation should provide to the study of self-regulated sys- 
tems, such as a unified description of the related phenomena, further comprehension of the internal 
structure and symmetries of the related equations, and the determination of the equilibria points via 
the energy function. We study the case of logistic systems, obtaining explicitly several s-equivalent 
Lagrangeans corresponding to the Verhulst's equation. Some dynamical properties are discussed in 
the light of this approach. Expressions for the mean energy function are also obtained. 
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I. INTRODUCTION 

The Lagrangean formalism is a very powerful tool in 
physics, presenting in a single formula all the dynami- 
cal information of a system. Nevertheless, it is widely 
accepted that a variational principle cannot be con- 
structed for an arbitrary differential equation Q : there 
exists a strict mathematical theorem that shows its 
existence for a given situation, and whose appHcation 
reduces the number of equations in physics that have a 
Lagrangean-Hamiltonian formulation ]^,|^. Due to its 
importance, efforts have been directed to make several 
methods to circumvent this condition and to construct 
modified variational principles Q, IE IB| - 

However, there is an interesting different approach, 
known as the inverse problem of the variational cal- 
culus, or IPVC (see, for example, 0, Q), that has 
its historical origin in the times of Helmholtz 0. It 
consists in studying the existence and uniqueness (or 
multiplicity) of Lagrangeans for systems of differential 
equations, meaning finding the Lagrangean, if it exists, 
from the equations of motion, instead of the traditional 
approach. This proceedure has become quite useful 
and one important result is that of Hojman et al. 10], 
who have proved that it is possible to construct the La- 
grangean for any regular mechanical system as linear 
combination of their own equations of motion. This 
particular construction, for example, is much wider 
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than the traditional definition L = T — V, which 
is only true when the "forces" involved are derivable 
from position-dependent potentials (or very few cases 
of velocity-dependent potentials) , therefore it may be 
used for general non-conservative systems. 

Whereas the equations of the models employed to 
describe auto-regulated phenomena can be understood 
as equations of motion in the variational sense, the 
latter approach shows itself to be one of the most ad- 
equate to obtain the Lagrangean formulation of the 
problem. Its appHcation to the study of self-Hmiting 
processes may provide additional understanding into 
the internal structure of these phenomena, also en- 
ables the use of a well known mathematical machinery 
to find conserved quantities, equihbria and stability 
cases, and other dynamical properties. 

In particular, the method proposed by Hojman et 
al. provides the desired Lagrangean both elegantly 
and without many complications, relating directly the 
constants of motion with the problem. Bearing this in 
mind, in the following pages we treat the problem of 
self-regulated systems described by the Verhulst's Lo- 
gistic Equation (VLE, from now on) [11, 12] by means 
of the Hojman proceedure, showing how a variational 
formulation can be easily obtained once we know the 
equation of motion and some additional information 
about the system. We have chosen the VLE because 
it had been well studied, thus we can easily compare 
the results here obtained to evaluate the approach, and 
also becasuse the model is simple enough as to give an 
idea, without involving too much time, of the way a 
self-regulatory system works. 

Historically this equation, firstly proposed by P.P. 
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Verhulst in 1838 as a continuous non-linear pop- 
ulation growth model with a self-limiting density de- 
pendent mechanism, was suggested as a direct gener- 
alization of the Euler and Malthus studies about 
population dynamics. 

If n(t) is the population at a time t, the rep- 
resentative time scale of response of the model to any 
change in the population and B define the carrying ca- 
pacity of the enviroment (see below), then the logistic 
equation can be written as 
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and we can inmediatly show that limt^oo n{t) = B (for 
all no), meaning, of course, that the carrying capac- 
ity corresponds to the the maximum size of the stable 
steady state population. In the Figure 1 the qualita- 
tive behavior of l(2} is presented. For no < B the 
profile is the characteristic sigmoid of the model, and 
for no > B the behavior is similar to an exponentially 
decay function. 

Indeed, Verhulst knew 0| that there were "biologi- 
cal obstacles" in the natural population growth, but to 
write a general mathematical expression for them was 
hardly possible, because each obstacle could affect in a 
different way the growth itself. Instead, he proposed a 
general "retardatrice" density dependent (cx n^) func- 
tion, a kind of compensating effect of overcrowding. 
However, it is important to remark that, even a very 
simple model, it is really useful as it matches a wide 
variety of different contexts, not only in population 
dynamics T^l, but also in biological growth (see for 
example, ^14]), radioactive decay series and even 
in some reaction-diffusion models fl^ . Effectively, all 
these cases (and more) can be adequately described by 
means of a characteristic saturation curve like ^ and 
in these cases an interpretation of the phenomenona 
by means of change in the population of some element 
(individuals, cells, isotopes, ...) can be made. 

The self-regulated systems, not only those described 
by the VLE, are present everywhere: certain aspects 
in stability of structures (firstly pointed out by Eu- 
ler ^3)) Townsend electron breakdown histeresis 
and magnetization processes the Amdhal law for 
scalability of computer programs and a long etcetera. 
As far as these phenomena are studied sometimes by 
quite different disciplines of science, it is interesting to 
explore if it is possible to provide a unified description 



for them in terms of families of Lagrangeans (or Hamil- 
tonians), that also would help to classify the systems 
by their dynamical properties. 

The distribution of this work is the following: Next 
section contains an abstract of the Hojman et al. 
method, which we use in Section 3 to construct the 
Verhulst's Lagrangian and Hamiltonian. Section 4 is 
devoted to the study of the energies related to the sys- 
tem, its equilibrium points and its constant of motion. 
Finally the concluding remarks are presented in Sec- 
tion 5. 



II. THE HOJMAN ET AL. METHOD 

For the purposes of this paper it is useful to work 
with the usually denominated second order approach 
In this context, the equation of motion for a 
mechanical system arises from a set of m differential 
equations: 

G' = q-r[q\q\t^^O, i,j = l,...,m (3) 

where the are the generalized coordinates and the 
point means total temporal derivative. 

In the IPVC the Lagrangean L ^g"' , , tj is con- 
structed such that relations (EJ can be effectively de- 
duced via the Euler-Lagrange equations. The exis- 
tence of such a Lagrangean is studied b^ means of the 
nowdays called Helmholtz conditions 
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Nevertheless, these conditions do not give any war- 
ranty about uniqueness. Two Lagrangeans are said to 
be solution-equivalents (or s-equivalents) if they just 
differ by a global multiplicative constant, 77, and a to- 
tal time derivative of some gauge A (^q'' , q^ , tj : 



dA 



(7) 
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The different systems of equations they provide, 
however, have exactly the same equations of motion. 
We can proceed now to describe briefly the method. 

It enables us to write L as a linear combination of the 
known equations of motion; then for i,j = l...m (the 
degrees of freedom) , 



L = 



(8) 



where 

Mi 



dq 
dA 

dq 
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and the are constants of motion of the mechanical 
system and the D2m-i are arbitrary functions whose 
arguments are constants of motion. One possible form 
for the D2m-i functions, given the conserved 
quantities, is presented in reference 0|. 

When the conserved quantities are unknown then 
the problem is reduced to find fXi such that the follow- 
ing system is satisfied: 



= 



(11) 
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where the on-shell derivative ^ is defined as 



dt dq' 



d 



d_ 

dt 



(13) 



For further details the reader is exhorted to review 
0, and the references therein. In the following 
section we make use of this proceedure to obtain the 
Lagrangean corresponding to the VLE. 



III. CONSTRUCTING THE VERHULST'S 
LAGRANGEAN 

For convenience, we shall write (QJ as ITU 



q=^ = kq{B - q) 



(14) 



with k = A/B, and we shall consider the one- 
dimensional problem. 

In our case, equation (EJ is 



q - kq{B - 2q) = 



(15) 



where q is given by II14|I . Explicitly in terms of the 
generalized coordinate (on shell), we have 

q-k^q{B - q){B -2q) = Q (16) 

together with the initial conditions 

g(0) = fcgo {B - qo) (17) 

9(0) = go (18) 

Note at this point that, indeed, the Verhulst La- 
grangean exist, for equation itTfill satisfies trivially the 
Helmoltz conditions. 

The second order Hojman et al. method provides 
then 



Lv=^i[q-k^q{B-q){B-2q)] 



(19) 



We just need to determine the factor ^, which in 
this case can be written as 



(20) 



There are several ways to proceed now. For exam- 
ple, observe that from the initial conditions and lll(i|l 
we can obtain the only constant of motion, C2, the 
system possesses (of course, for m degrees of freedom 
there are (2m — 1) functionally independent constans 
of motion |l21|); the arbitrary function C\ has then a 
prescribed form Q, and fj, is obtained by means of 
ll271|l . Furthermore, it is possible to solve lfTT|l in the 
one-dimensional case for fi, which yields a general so- 
lution in terms of Bessel functions. Using the initial 
conditions of the problem we find the particular solu- 
tion. It is also possible to introduce the expression for 
the Lagrangean in the Euler-Lagrange equation and 
compare with the equation of motion to obtain /J.j2fl|- 
We shall ilustrate here the first way. 

It is not difficult to show that 



C2 ^ -g - -k'q' {B -qy^O (21) 
and a possible choice for Ci is (see reference 01) 



(22) 



4 



where c is an arbitrary constant that multipHes the 
equation of motion. The Lagrangean reads then 

Ly = cClq [q - k^q{B ~ q){B ~ 2q)\ (23) 

When inserted into the Euler-Lagrange equation, 
^2'^\ provides the corresponding equation of motion 
(I16|l plus terms that are zero in virtue of the constant 
of motion. However, it has been shown that for the 
one-dimensional problem there exists an infinite num- 
ber of Lagrangeans Thus, we can write a more 
simpler s-equivalent one by means of the total time 
derivative of a certain gauge. Given the constant of 
motion, let choose 



rfA _ .2 /„ dC2 



(24) 



Without loss of generality we set c = 77 = 1, and 
together with Ij24|l . equations JJjl and i2'M provide 



Lv^lq' + hViB~qf 
and the Hamiltonian can be written as 



Hv = lp'-lkV{B-qf 



where the generalized momentum is 

dq ^ 



(25) 



(26) 
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IV. CONSERVATION, ENERGIES AND 
EQUILIBRIUM CONDITIONS 

Now we proceed to study some of the applications 
the Lagrangean formalism offers to the Verhulst sys- 
tem. 

The first term in II25|I and Il2(ill is identified as the 
usual traslational "kinetic" energy: quadratic and ho- 
mogeneous in the first temporal derivative of the gen- 
eralized coordinate. The second term on the right side 
of equation ()2(ill is explicitily time- independent, mean- 
ing that Hv coincides with the energy function and it 
is then possible to write the potential energy of the 
system as 



V 



■hWB-qf 



(28) 



The behavior of this potential is qualitatively pre- 
sented in Figure 2 for fixed carrying capacity. In 



this potential, the smaller the time scale (1/^), the 
more closed the curve and the larger the depression. 
Likewise, the greater the carrying capacity, i?, the 
more pronounced the depression and the more open 
the curve. Equation II28II can be interpreted as an 
impulsor-retardatrice potential, depending only on the 
generaHzed coordinate. In fact, observe that when 
qa < B, for the interval (0, B/2) the acceleration asso- 
ciated impulses the movement, being conversely for the 
interval {B/2,B); the point B/2 corresponds to a lo- 
cal minimum, where the acceleration instantaneously 
annuls itself. On the other hand, when qo > B, i.e. 
the interval {B,oo), the character of the acceleration 
is always impulsive, even being g < 0. All this is con- 
sistent with the profiles sketched in Figure 1 and with 
the standard knowledge about equation ifnji . 

Despite the fact the VLE describes multiple phe- 
nomena which one should think presents dissipation of 
energy (for example as in population dynamics, growth 
of living beings or meteorological models), it is easy to 
show that the Verhulst system is conservative (popula- 
tion systems with a first integral has been studied, see 
[T^ and references therein); thus, the energy function 
Hv equals the total energy E. Note that in virtue of 
(I14II equation (I26II is identically zero. Taking all these 
considerations into account, it is clear in Figure 2 that 
the kinetic energy tends to annul itself as q approachs 
B. As a consequence, the system takes an infinite time 
to reach the steady state, and qo < B implies always 
q < B (similarly for go > B). 

Another interesting issue is that there are not oscil- 
latory behaviours, despite the similarity with harmonic 
or the anharmonic potentials. Observe that as far as 
we are just interested in feassible results, we must take 
only the positive values of q{t). Thus the character of 
the force (see equation is not restorative, but the 
commented before for each interval. 

One important application is that of finding the sta- 
bility cases. The usual derivative criteria allows us to 
find the equilibrium points, for which 



dE 



= 



The stability is then ruled by the condition 

d^E „ 
> 



(29) 



(30) 



The equilibrium solutions are qe — {O, The 
first one corresponds to a point of instability, the sec- 
ond to minimal potential energy and the last one is, 
as expected, the state of stationary equilibrium. As 
before, this analysis is in accordance with the known 
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behaviour of the VLE (see Figure 1). It is important 
to remark here that usually the study of its equilib- 
rium points is treated by means of a Taylor expansion 
of ifnt about q = and q = B (see for example [T^V 
The variational formalism provides thus both a more 
elegant and complete treatment. 

On the other hand, as far as the VLE can model 
various self-regulated phenomena, it should be useful 
for the description of the related statistical systems to 
write the total mean energy. We shall use the virial 
theorem for this porpouse. In our case. 



(T) 



1 /dV 



2\dq 



e {q^ (B -q){B- 2q)) 



(31) 



Then 



B 



Note that if 
to 



4g 



Bk + ^jB^k^ - 4fcgJ J + 2{E) (37) 
^ 1 the latter equation simplifies 



.2\ AB , , 

q-) = —{q)+2{E) 



(38) 



The use of the other solution, q-{t), provides 
formally 

(E) = (39) 
and then equation l(37jl is no longer true. 



Then, the total mean energy is written as 

{E)^ {T + V) =-h^q\B-q)) 



kid 



6\dt 



(32) 



or, explicitily, evaluating the integral between t — 
and t ^ oo, 







(33) 



Thus the total mean energy of the systems described 
by VLE is always null. But if as superior temporal 
bound we choose the characteristic time of the system, 
then 



6B 



4 



(34) 



Returning to the general case, observe that under 
certain circumstances equation l(32|l make it possible 
to stablish a relationship between the mean-squared 
velocity and the mean-linear velocity. Indeed, solving 
l|14ll for q we find two solutions: 



kB ± y^B^k^ - Akq 



2k 



Using q+{t), equation II32II can be restated as 
/ ■ 2 



(i;) = / y - j<? [Bk + ^jB^k^ - Akq 



(35) 



(36) 



V. CONCLUDING REMARKS 

We suggest in this work the use of the IPVC to ob- 
tain a variational formulation for self-limiting systems. 
This approach provides at least a Lagrangean, if it ex- 
ists, departing from the equation of motion and regard- 
less the system being conservative (when the definition 
L = T—V holds) or not. There are several advantages, 
namely, it is possible 

(a) to have in just one relation all the dynamical 
information of the system; 

(b) to use a well know mathematical machinery to 
obtain from the Lagrangean or Hamiltonian the con- 
served quantities (and symmetries), equilibria points 
and trajectories, etc.; 

(c) to group self-regulated phenomena having the 
same variational formulation, even when describing 
quite different scenarios, classifying the systems by the 
form of the Lagrangean or by its dynamical properties; 

(d) to solve, when possible, the equation of mo- 
tion. For some systems (as those described by Lotka- 
Volterra models) the Hamilton-Jacobi formalism, for 
example, should be useful to integrate the equation of 
motion; 

(e) to help to understand how the models are re- 
lated. Some models are, indeed, particular cases of 
more general equations; 

(f) to test if a certain model is acceptable in terms 
of the Least- Action Principle. 

On the other hand, making use of the Hojman 
et al. method, we have found the Lagrangean ll25|l 
and Hamiltonian l|26ll corresponding to the (1+1)- 
dimensional logistic equation as an example of the 
way a variational formalism can be obtained for gen- 
eral self-regulated systems. This may be regarded as a 
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proof that the VLE corresponds to an extremal trajec- 
tory, thus it is physically acceptable (it complies with 
the Hamilton's Least- Action Principle) despite the fact 
Verhulst introduced it heuristically. 

The potential energy II28|I was briefly analized and 
the equilibria points and the stable/unstable cases 
were obtained by means of the usual derivative criteria, 
in accord with the other method (see 0|). Also, using 
the virial theorem it was possible to write an expres- 
sion for the total mean energy of the system, equation 
1I32I1 , and a relation between the mean-squared velocity 
and the mean-linear velocity, equation ll37ll . 

It is not surprising to note that the Malthusian limit 
is trivially obtained for the particular case when B 
tends to 00. Thus, the potential energy for this case 
can be written as 

VMalthus = -^^'9' (40) 

We have in this work studied the thus-called sec- 
ond order approach. There is, certainly, a flrst order 
method of the Hojman et al. proceedure 0, E|; we 
show briefly in the Appendix the way an s-equivalent 
Lagrangean for the VLE can be constructed using this 
approach. Even though the Lagrangean reproduces 
the logistic equation, as required, the physical inter- 
pretation of its terms is not so clear. However, gen- 
erally speaking, it is important to bear in mind this 
proceedure because it sometimes provides a variational 
formalism of the physical phenomenon while the sec- 
ond order approach does not 

Finally, it could be useful to study in future works 
generalized auto-regulated models as for example 
those reported in references fH il El EE El]. It 
is also interesting to extend this work to the study of 
the quantum systems related with self-regulated phe- 
nomena l23l. 
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APPENDIX 

In order to obtain an s-equivalent Lagrangean for 
the Verhulst 's equation using the flrst order approach 
of the Hojman et al. method we proceed as follows. 



First, let write the logistic equation as 

q = kq{B -q) 



(41) 



where q = qit) is the generalized coordinate. Then, let 
write and choose a new integral variable Q{t) 



(42) 



Q{t) = J q{t)dt 
such that equation lUTll transforms in 

Q{t) = kQ{t) \b - Q{t) 



(43) 



This second order equation can be casted in two flrst 
order equations 



Xi = X2 
X2 = kx2{B - X2) 



(44) 



where Xi = Xi{t) and X2 = X2{t). Note that indeed X2 
is exactly the generalized coordinate q. 

The solution of this system of equations is given by 



xi ^ / X2dt + Ci 



B 



X2 = 



1 + BC2e-^Bt 



(45) 



(46) 



and this functional form for X2 is just another way to 
write the usual form JSJ. Now, solving and iTIfi)) 
for Ci and C2 we obtain 

Ci = kxi - kBt -InB -\nx2 (47) 

^kBt 



C2 = 



X2 



(48) 



and following Hojman and Urrutia the s-equivalent 
Lagrangean is flnally written as 



kBt 



kxi kxi 
~2 B~ 
kB\n{x2) , ,2 



kB 

X2 



X2_ 

■''2 



X2 



X2 



-\- k xi + fclnn(a;2) -|- 



BX2 

kx2Xi 



k^Bxi k'^BH kB\n{B) 

X2 X2 X2 

kBtx2 \n{B)x2 ln(a;2)i2 



(49) 
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Figure captions. capacity {B = 20). 

1. The Verhulst logistic equation for fixed carrying 2. The Verhulst potential for differents k {B 
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